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ABSTRACT 

Demand  during  lead  time  is  considered  as  the  sum  of  a  random  number 
of  random  variables.  Discrete  random  variables  are  considered  as  de- 
scribing demand  and  lead  time  and  probability  generating  functions  are 
used  to  determine  the  probability  generating  function  of  demand  during 
lead  time.  The  probability  distribution  of  demand  during  lead  time  is 
extracted  from  the  probability  generating  function  by  using  a  power 
series  expansion  and  the  calculation  of  probabilities  from  the  expan- 
are  discussed.  The  mean  and  variance  of  demand  during  lead  time  ar© 
also  obtained  from  the  probability  generating  function,. 
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1.   Introduction 

The  two  major  problems  in  any  inventory  control  system  are  to 
decide  upon  the  amount  to  order  each  time  replenishment  is  necessary 
and  to  decide  what  should  be  the  minimum  stock  level,  of  any  items 
at  which  an  order  should  be  placed.  The  second  problem  is  the  subject 
of  this  paper. 

The  minimum  stock  level  at  which  an  order  for  replenishment  of  aia 
item  is  the  reorder  point,  denoted  by  R.  The  reorder  point  is  a  func~ 
tion  of  the  demand  and  the  lead  time  of  an  itemu  The  demand  for  an 
item  is  the  rate  at  which  the  stock  of  the  item  is  depleted.  The  lead 
time  of  an  item  is  the  length  of  the  interval  from  the  time  a  replenish 
ment  order  is  placed  until  the  order  is  in  stock,  ready  to  satisfy  the 
customer's  demand, 

A  simple  deterministic  model  of  an  inventory  system  with  constant 
demand  and  lead  time  is  shown  in  figure  1, 
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Deterministic  Inventory  Model 
Figure  1 


The  stock  is  depleted  at  a  uniform  rate  from  a  stock  level  of  S 
until  it  reaches  the  reorder  point.  When  the  stock  level  of  R  units 
is  reached,  a  replenishment  order  for  S  units  is  released  and  the 
order  is  received  at  time  A,  just  as  the  stock  is  exhausted 0  The 
stock  replenishment  raises  the  stock  level  back  to  S  units  and  the 
process  is  repeated. 

Rarely  is  demand  uniform,  however,  and  there  is  uncertainty  about 
when  the  R  units  of  stock  on  hand  at  the  time  of  placing  an  order  will 
be  exhausted.  A  more  realistic  model  of  an  inventory  system  with 
variable  demand  and  constant  lead  time  is  shown  in  figure  2„ 
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Variable  Demand  Inventory  Model 


Figure  2 


S,  R  and  A  are  the  same  in  figure  2   as  they  were  in  figure  1. 
The  expected  depletion  of  stock  is  along  the  line  S  A.  If  the  lead 
time  was  deterministic,  an  order  placed  at  R,  time  X,  would  be  re- 
ceived at  time  A,  In  the  time  interval  during  which  replacement 
stock  is  an  order,  X  A,  there  can  be  a  variation  in  the  demand j 


therefore,  the  demand  may  deplete  the  R  units  of  stock  along  any  line 
between  a  maximum  demand  >  R  A1,  and  a  minimum  demand  5,  R  A",  shown  by 
the  shaded  area  A'  R  A"0  The  variation  in  demand  allows  for  the  pos~ 
sibility  of  depleting  the  stock  at  time  A1  and  not  being  able  to  meet 
demands  for  stock  in  the  period  A'  A.  There  is  also  the  possibility 
that  the  demand  could  be  as  low  as  R  A"  and  in  this  case  there  would 
be  B  units  of  stock  on  hand  when  the  replacement  order  is  received,, 

A  third  and  the  most  realistic  model9  considering  both  variable 
demand  and  lead  time,  is  shown  in  figure  3.  The  first  two  mode  La 
were  special  cases  of  this  model  and  therefore  model  3  will  serve  to 
describe  the  inventory  system  considered  throughout  the  rest  of  this 
paper. 
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Variable  Demand  and  Lead  Time  Model 


Figure  3 


The  variation  in  the  lead  time  will  cause  the  replenishment 
time  interval  X  Y  to  vary.  Time  X  is  the  time  at  which  the  reorder 


point,  R,  is  reached,  but  it  is  no  longer  certain  that  the  order 
will  be  received  at  time  Y„  The  order  may  be  received  at  time  U  or 
time  V  or  any  time  in  between,  if  we  let  time  U  and  time  V  represent 
the  minimum  and  maximum  times  of  receipt  of  the  order 0  Figure  3  is 
similar  to  figure  2,  but  the  zero  stock  level  of  figure  2  is  the 
same  number  of  units  below  R  as  C  is  in  figure  3.  Figure  3  then  is 
the  same  diagram  except  that  the  entire  stock  is  increased  by  C  unit.; 
to  allow  for  the  variation  in  lead  time  as  well  as  the  variation  in 
demand.  Let  us  consider  the  four  possible  extreme  cases  of  maximum 
and  minimum  demand  and  lead  time8  A  replenishment  order  is  placed 
when  the  stock  on  hand  is  R  units  at  time  X„  If  the  maximum  demand 
and  the  minimum  lead  time  occur,  the  stock  will  be  depleted  along 
line  R  A1  until  there  are  D  units  at  time  U  when  the  replacement  stock 
is  received.  If  the  minimum  demand  and  the  maximum  lead  time  occur ^ 
the  stock  will  be  depleted  along  line  R  A"  until  there  are  E  units  in 
stock  at  time  V  when  the  order  is  received e  Minimum  demand  and  mini- 
mum  lead  time  will  result  in  F  units  of  stock  on  hand  at  time  U  when 
the  order  is  received.  In  each  of  these  three  cases 9   R  units  of  stock 
have  been  more  than  enough  to  provide  protection  against  stockouts 
with  both  variable  demand  and  lead  time.  Only  in  the  case  of  maximum 
demand  and  minimum  lead  time  would  the  stock  be  below  C  units „  The 
last  case  is  the  combination  of  maximum  demand  and  maximum  lead  time0 
The  stock  will  be  depleted  along  the  maximum  demand  line  R  A"  and  will 
be  exhausted  at  time  Z.  The  order  will  not  be  received  until  the  maxi- 
mum lead  time,  time  V,  therefore  the  system  will  be  out  of  stock  for  a 
period  of  V-Z  time  units. 


From  figure  39  it  can  be  seen  that  if  the  demand  follows  the 
line  R  A  and  time  Y  is  the  lead  time>  then  the  demand  during  lead 
time  is  equal  to  R  «  C  units  of  stock,,  It  can  also  be  seen  in 
figure  3  that  as  the  height  0  C  is  increased  the  intersection  of 
the  actual  demand  and  the  actual  lead  times  whatever  they  may  be,, 
will  intersect  above  the  zero  stock  level  more  often  than  below. 
The  reorder  point,  R,  is  then  the  sum  of  the  expected  demand  daring 
lead  time  plus  a  safety  level  or  buffer  stock  of  C  units 0  The  re- 
order point  equation  is  then; 

R  =  Z  +  V  S  L  (1.1) 

where 

R  =  reorder  point 

Z  =  expected  demand  during  lead  time 

V  S  L  =  variable  safety  level 

The  safety  level  is  denoted  as  a  variable  safety  level  because  the 
protection  against  a  stockout  increases  as  the  number  G  is  increased,, 
The  safety  level  can  be  varied  so  that  no  stockouts  occur  or  soma 
small  number  of  stockouts  occur  in  a  specified  number  of  replenish- 
ment periods ,  That  is,  the  inventory  manager  may  decide  that  a 
safety  stock  large  enough  to  protect  against  no  stockouts  would  be 
too  expensive  and  therefore  he  is  willing  only  to  maintain  a  safety 
stock  that  will  protect  against  stockouts  ninety  or  ninety-five  per=> 
cent  of  the  time.  The  safety  level  is  also  a  function  of  the  reorder 
quantity  since  the  larger  the  reorder  quantity  for  a  given  lead  time 
the  less  frequently  the  system  is  exposed  to  a  stockout.  This 


consideration  is  taken  into  account  when  the  protection  level  is 
determined  and  does  not  affect  the  reorder  point  equation  J6j  9  111 


2.   Discussion  of  the  Problem 

The  availability  of  alternative  courses  of  action  implies  the 
need  for  managers  to  make  decisions 0  The  selection  of  the  reorder 
point  for  an  item  in  the  manager's  stock  is  then  a  decision  problem 
The  reorder  point  rule  must  satisfy  an  objective  function  (such  as 
maximum  service)  subject  to  some  constraints  imposed  by  higher 
authority,  available  funds ,  available  facilities 9   etc0  In  almost 
all  cases,  the  method  by  which  the  problem  is  solved  requires  a  model 
of  the  operations „  Such  a  model  is  based  on  a  set  of  relations  link- 
ing the  important  variables w  Figure  3  is  the  geometric  model  of  tha 
inventory  system  being  considered  and  the  reorder  point  equation  (Id) 
will  lead  to  the  mathematical  model  which  will  be  used  to  determine 
the  probability  distribution  function  of  demand  during  lead  time* 

The  reorder  point  equation  is  a  rule  by  which  the  inventory 
manager  may  decide  when  to  reorder  replacement  stocku  The  variable 
Z  is  the  mean  of  the  probability  distribution  of  the  demand  during 
lead  time  and  the  variable  V  S  L  is  a  function  of  the  variance  of  th© 
same  probability  distribution,,  Therefore  $  in  order  to  use  the  reorder 
point  equation,  the  probability  distribution  function  of  demand  during 
lead  time  or  at  least  the  mean  and  variance  of  the  distribution  func- 
tion must  be  known , 

Let  Z  be  a  random  variable  which  describes  the  demand  during 
lead  time j,  then  the  number  of  demands  that  occur  during  a  lead  time 
may  be  written  ass 

Z  =  D1  +  D2  +  909  +  DL         (2d) 


where 

Z  =  the  number  of  items  demanded  daring  a  lead  time 
D.-  demands     ,    i  =  1,  o»o  9  L 
L  =  lead  time 

It  is  assumed  that  the  demands  are  independent  and  identically  dis- 
tributed and  demand  and  lead  time  are  independent  discrete  random 
variables  with  known  probability  mass  functions 0  In  general s  this 
paper  will  only  consider  discrete  random  variables  since  demands 
naturally  occur  at  discrete  intervals  and  lead  times  are  normally 
described  in  days,,  weeks  or  months0  Equation  (201)  is  then  a  model 
which  describes  the  demand  during  lead  time  by  linking  together  the 
random  variables  D  and  L  and  is  a  model  that  can  be  used  to  deter- 
mine the  expected  demand  during  lead  tirne^  Zs   and  the  variable 
safety  level y  V  S  Ls  Equation  (2d)  can  be  described  more  concise- 


ly as 


2  =  U  Di  (2-2) 


Using  equation  (2e2)  as  the  model  of  demand  during  lead  time,, 
it  will  be  shown  that  the  mean  and  variance  can  be  computed  for  any 
probability  mass  functions  that  describe  demand  and  lead  time  or 
from  the  mean  and  variance  of  demand  and  lead  time  if  that  is  all 
that  is  known  about  them*  It  will  also  be  shown  that  the  distri- 
bution of  demand  during  lead  time  can  be  found 5  in  the  discrete 
cases  by  the  use  of  probability  generating  functions 0  After  the 
probability  distribution  is  found s  it  will  then  be  shown  that  it 
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can  be  used  to  set  the  reorder  point  using  a  variable  safety  level j 
given  the  level  of  customer  service  desired  0  The  final  section 
treats  a  special  case  in  which  lead  time  is  a  continuous  random 
variable  and  the  resulting  distribution  of  demand  during  lead  time 
has  a  familiar  negative  binomial  form. 


3»   The  Mean  and  Variance  of  Demand  During  Lead  Time 

The  mean  and  variance  of  Z  are  derived  using  general  probability 

mass  functions  for  D  and  L0 
Lets 
p  (d)  be  the  mass  function  of  demand  with  mean  D  and  variance 

p  (1)  be  the  mass  function  of  lead  time  with  mean  L  and  vari- 
ance U 

p  (z)  be  the  mass  function  of  demand  during  lead  time  with 
mean  Z  and  variance  CJ_ 

The  mean  of  Z  is  the  expected  value  of  Z  and  is  denoted  as  Z 

E(Z)  =  E  [e(zJL)] 
=  E  [l  E(D) 
=  E(D)  E(L) 
=  D  L 

Therefore  the  expected  demand  during  lead  time  Is 

Z  =  D  L 

The  variance  of  demand  during  lead  time  Is  derived  as  follows  § 


Var(Z)  =  E  [var(Z|L)J      +  Var     |e(Z|L)1 
Var(Z)  =  E   [_L   6^  J     +  Var     |_L  E(D)J 
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E(L)  Op2  +   [E(D)j  <Jl 


l  a2  +  d2  cr2 

D         L 


The  mean  and  variance  of  demand  during  lead  time  are  functions 
of  the  means  and  variances  of  the  demand  and  lead  time  only,  there- 
fore a  knowledge  of  the  mass  functions  of  demand  and  lead  time  is 
sufficient  to  calculate  the  mean  and  variance  of  demand  during  lead 
time.  They  may  also  be  calculated  empirically  using  historical  de= 
mand  and  lead  time  data  to  estimate  the  means  and  variances  of  demand 
and  lead  time0  The  mean  and  variance  of  a  distribution  function  may 
also  be  calculated  from  the  probability  generating  function  of  the 
distribution.  The  probability  generating  function  of  a  probability 
distribution  function  is  discussed  in  the  next  section  and  the  calcu= 
lation  of  the  mean  and  variance  from  a  probability  generating  func- 
tion is  shown  in  Appendix  I„ 
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4c   Distribution  of  Demand  During  Lead  Time 

In  this  section,  the  probability  generating  function  of  demand 
during  lead  time  will  be  derived  in  general „  The  generating  func- 
tion will  then  be  used  to  obtain  the  probability  distribution  of 
lead  time  demand  by  expanding  the  generating  function  in  a  power 
series  and  extracting  the  coefficients 0  Probability  generating 
functions  are  defined  for  discrete  random  variables  only0  There- 
fore j  the  discussion  will  be  limited  to  discrete  random  variables, 

The  probability  generating  function  of  demand  during  lead  time 
is  the  generating  function  of  the  sum  of  a  random  number  of  random 
variables  and  is  defined  as 


y  k  «  r  ,i 


h(s)  =  E(sZ) 


To  derive  h(s),  let  the  probability  generating  function  of  demand 
and  lead  time  be 


g(s)  =  E(s^) 


f(s)  =  E(sD) 

and 

.Z\  _  t?  i  JL  It  - 


now 


h(s)  =  E(s*)  =  E  (s*  L  =  1) 


E(sZ|L=l)  =E(sD^   sDi  oousDl) 


=  E(8D|)E(sPt)...E(sDl: 
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since  the  Dj_  were  assumed  independent  and  identically  distributed 
Also  since  f(s)  =  E(s  ),  we  have 


and 


since 


E(sZ|L=  1)  =  [f(s)]1 
E(sZ|L)  =  [f(s)JL 

E  [e(sZ|l)]  ==E  j[f(s)]] 
=  g  [f(s)] 
g(s)  =  E(sL) 
,\  h(s)  =  g  [f(s)] 


Thus,  if  the  distribution  functions  of  demand  and  lead  time  are  known 
or  assumed,  the  probability  generating  functions  of  demand  during  lead 
time  can  be  derived.  Then  with  a  knowledge  of  h(s),  the  probability 
distribution  of  demand  during  lead  time  is  known.  Examples  showing 
the  computation  of  the  distribution  of  Z  will  be  given  later , 

Implementation  of  the  model  given  in  equation  (2.2)  may  generally 
take  place  using  familiar  discrete  distributions  for  the  variables , 
since  such  distributions  have  convenient  mathematical  properties. 
These  familiar  distributions  usually  have  positive  probability  at  zero 
and  this  situation  usually  is  not  compatible  with  lead  time.  Prac- 
tically, lead  time  can  never  be  zero,  if  for  no  other  reason  than  it 
takes  some  time  to  prepare  an  order.  On  the  other  hand,  if  the  data 
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in  the  inventory  control  system  is  based  on  the  month  as  the  basic 
time  units  a  lead  time  of  one  day  could  be  considered  zero.  In  any 
event,  if  the  mean  of  the  lead  time  distribution  is  small,  then  "the 
use  of  a  popular  distribution  will  result  in  probability  at  zero  which 
is  not  negligible.  Mathematically ,  the  model  is  correct  since*  when  L 
equals  zero  the  sum  is  empty  and  Z  is  zero.  Demand  during  lead  time 
can  also  be  zero  when  demand  is  zero  during  the  lead  time  or  when 
both  demand  and  lead  time  are  zero. 

There  are  several  methods  by  which  the  input  distribution  can  b© 
adapted  to  consider  either  allowing  or  not  allowing  some  probability 
for  lead  time  to  be  zero.  In  the  first  ease,  when  lead  time  can  be 
zero,  the  straight  forward  use  of  the  model  with  distribution  func- 
tions of  demand  and  lead  time  that  are  defined  on  the  non-negative 
integers  will  have  provided  the  desired  results a  There  are  two  methods 
that  may  be  used  to  prevent  the  consideration  of  zero  lead  time  in  the 
model,,  These  methods  redistribute  the  probability  mass  functions  from 
a  range  over  the  non=negative  integers  t©  a  range  over  the  positive 
integers.  The  first  is  translation,  i.e.,  let  L  =  1+X,  where  X  takes 
non-negative  integral  values?  and  the  second  is  truncation,  i  ec, 
use  a  frequency  function  that  is  conditioned  on  the  variable  being 
positive.  Let  us  consider  the  Poisson  distribution  as  an  example 0 
The  Poisson  distribution  is  usually  defined  ass 


p(xj  A)  =  2 — A.        for  A>0,  x 
x! 


The  mass  function  can  be  translated  over  the  positive  integers 


H 


writing  it  as: 


-A  oc-1 

p(x;  X)   -  - 4 ,   f or   A  >  0,  x  =  1,2, 

(x-1)! 


or  the  distribution  function  can  be  truncated  by  using  a  conditional 
probability.     That  iss 


P  [x  =  x|x  >  o]   =  e 


x!(l-e~A) 


1  A 


eA  -1  x! 


Ax 

p(x;M  =  ,  for  A  >  0,  x  =  1,  2,oo 

(eA-l)  x! 


There  are  then  three  choices  offered  for  adjusting  the  input  to  fit 
the  situation  or  the  preference  of  the  user. 

Let  us  now  consider  some  generating  functions  of  demand  and  lead 
time  so  that  they  will  be  available  for  use  in  determining  the  prob- 
ability generating  function  of  Z.  The  Poisson  and  negative  binomial 
distributions  have  frequently  served  as  good  approximations  for  the 
distribution  of  demand  so  they  will  be  used  in  the  following  examples 0 
Lead  time  distributions  have  not  been  as  well  documented 9   but  the 
negative  binomial  distribution  is  a  good  model  for  discrete  waiting 
time  random  variables  and  when  the  parameter  r  equals  one,  the  nega<= 
tive  binomial  becomes  the  geometric  distribution.  As  a  second  lead 
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time  distribution j  the  Poisson  distribution  will  be  used0  The  use  of 
the  Poisson  distribution  for  lead  time  may  not  be  a  good  example 9 
with  a  large  mean,  the  Poisson  distribution  roughly  approximates  & 
normal  distribution  with  the  mean  equal  to  the  variance  f  l"]0 

The  generating  functions  for  the  Poisson  and  negative  binomial 
distributions  in  each  of  the  three  cases  are  given  below „ 
For  the  Poisson,  when 

p(x$  A  )     =    —  ■  ■'"  s    X    =    0919     oco 

The  probability  generating  function  iss 

6Q> 


f (s)  =  E(sX)  =  Zj    sk  p[x  =  kl 


=  e 


k=0 
a(s-1) 


Ud) 


when 


e~*     XX=1 
p(x|  A)  =  ~. —  -Q- —  j,  x  =  1,  2, 


(x-1)! 


o  e  a 


,    N  A(s~*l) 

f(s)  =  s  e  (^02; 


when 


X™-       I  jj      <£  ^        ooo 


for  negative  binomial  distribution,  whes 
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n(w.  „\  _  fr+x-1\     r  nx 
pU$r,p)   =  I  I   p     q 


^X~Uft  I     y  CO© 


f  (s)  =  f— Mr  U.3) 

\1-q  s  / 


when  the  distribution  is  translated 


p(x;r,p)  =  (r+*-Z)   pr  q*"1  ,1  =  1,2,.,, 


f(s)  =  ■  P    \F 


= •  (— y 

\1-qs     / 


when  the  distribution  is  truncated 

(r+x-f 


r     x 

p(x$r,p)  =   V    X    '    v  ,  x  =  1,  2, 

1-P 


f(.)  =  pr  Qjig£jl 

i-Pr 


As  an  example  of  determining  the  probability  generating  fune= 
tion  of  demand  during  lead  time,   let 

D  be  a  Poisson  random  variable 

L  be  a  translated  negative  binomial 
then 

f(s)  =e  Ms~1) 
and 

h(s)  =  g  [f     (s)] 
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<a& 


I 


r     i_    e  A(s-1)k     /r+k~2\      k  1 
k=1  I     k=li 


-  ^r  .  A(i 


=    p       8 


Pre 


=1)  z. 


-(s-1)    (k-1)  /V+k~2\      k-1 


k=1 


k-r 


A(s-D     2      (r+M    L  / 
k?i      \     k-V    L 


*(s~l) 


k-1 


pr  e  A(e-1) 


1-qe 


The  probability  generating  functions  of  demand  during  lead  times   Z9 
were  derived  since  they  can  be  used  to  find  the  probability  of  de- 
mand during  lead  timee     That  iss   since 


-I 


h(s)   =  E(s*)  =  Z-i    s"  P   |Z  =  k 

k=0 


Then  in  the  power  series  expansion  of  h(s)9  the  coefficient  of  s 
is  the  P  Z  =  k I  .  For  example,  writing  out  several  terms  of  h(s)9 
we  haves 

h(s)  =p[z  =  o]+spE-i|    +s2pfz  =  2J    +eoo 


It  therefore  is  necessary  to  put  h(s)  in  the  form  of  a  power  series 
in  order  to  determine  the  probability  of  demand  during  lead  time 
for  any  k. 

Using  the  probability  generating  functions  for  the  variable  Z 
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in  the  previous  examples,  we  can  put  them  in  a  more  useful  form, 
For  Poisson  demand  and  translated  negative  binomial  lead  time 


h(s)  =  e   X(s-1) 


JL 


__1-qe 


M 


»1)J 


OO 

=eA(s-1)prS(rr)[-A(s-i)J 


oo 


=  pr 


=  pr 


Zj     (r+k-1^ 

k=o  V   k   / 


qk  esMk+l)  -  (k+1) 


oo  co 

2  fr+k-A  qk  .-A(fcfi)  2  s"  J^teafl! 

k=0  \     k    /  n=o  n! 

OP  oo 

=  IsnprI    (rfk-l\    k  e-  A(k+1)  j>lk+l2P 
n=0  k=0    V     k   /  n! 


and 


.  P     Z  =  nj  is  the  coefficient  of  sn  in  the  power  series 


CO 


rz=n]=Pr2(rr)^«-A(k+i)lAiilI]1! 

L  J  k=0  V     K    /  n! 


-  p    e 


oo 

4^2    (r+H    (qeA)k   (k+On 
n!      k=o  \     k    ; 


The  calculation  of  these  probabilities  can  be  reduced  if  the  sum 
from  zero  to  infinity  is  reduced  to  a  finite  number  of  terms  with  a 
very  small  remainder.  A  method  for  finding  the  number  of  terms  to 
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be  summed  until  the  remainder  is   negligible  will  be  shown  along  with 
a  computing  technique  for  determining  the  probabilities. 
To  find  a  bound  for  the  remainder  of  the  series,   let 


and 


then 


AK 


00  r— i 

=  Zuu  =  Z>+k-l\     (     -a  )k  (k+l)n 
k=0  k       k=0\    k     j 


a  -  qe 


CO  CO 

"K+1+1 


7  7  Uf 

^       ulr   =    %+1     ^    — T 


RlC  "  Hf1   "k  ""  "K+1 


UK+1+1 


-  (Shi)  ak+1+l  (k+2+Dn 

Uk+1  ££)     aK+1      (k+2)n 

(r+K+l)l(K+l)l         /K±2±l\n  a1 

(K+1+1) !(r+K)!         I    K+2  / 


(r+K+1)    (r+K+1-1)...    (K+l+r-(r-2)    )     /K+2+l\n         1 
(K+r)    (K+r-1)...      (K+r-(r-2)    )  V  K+2  1 


=     (r+k+1)    (r+K+1-1)...    (K+l+2)         /K+l+2\n       &1 
(K+r)(K+r-l)    ...    (K+2)  V  K+2  / 


Dividing  term  by  term 
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1 

+  — 
K+r 


)  (i+_U  ...  (1+-LV(1+-L\V 

/    V         K+r~1/  V         K+2/   \         K+2/ 


(•4-) 


a1 


then  co 

n+1 


^^"k+i  1=6  (1  +t)       al 


It  is  possible  to  dominate  the  right  hand  side  with  a  integral  pro= 
vided  the  integrand  is  monotone  decreasing.  This  will  be  true  if  K 

is  sufficiently  large „  Thus,  let 


=>B  =  la  a  -  -  A  +  In  q 
and 

¥(D « (i+-jf) 


l\n+1      -Bl 


then 


■  <!>t)" •-'  [- ("t) ■ ' f ] 


Setting  the  derivative  equal  to  zero,  we  have 


1  =  £±1-K 


n+1 

. * „  the  maximum  for  1  occurs  if  K  >  — — 

B 

hence  we  must  use  K  such  that  it  is  at  least  this  large . 

n+1 


far  K  > 


B 
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iJO  CO 

&6   V  K  /  fci     V  K  / 


n+1   , 
^1+1        (1  +  ■— J       e         dx 


I  (T)  ©'  I  ■>  •- 


e"5*  dx 


making  a  change  of  variables,  let  v=Bx,  then  dv=Bdx,  and 

<s©  n+1  poo 

z  (. .  4-r  •• « •  •  i  (-1  ©'  *  i .« 


e       dv 


3=0 
n+1 


I ffl  S1  ^ ' 


B    1     + 

3 


(n+1)!      /  1\n+1 

^1+       B 


(1  +  -IV 

V         KB/ 


which  is  a  bounded  function  of  K,  thus 


"kc"k+i    L1  + 


a.-<-w  [i+^  (i^)n+1] 


now  if       KB  >n+1 
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%<(£*)     e=B(K+D    (K+2)n     [l  +  <S±Hi     .] 

In  practise 9  all  of  the  constants  and  parameters  in  the  above 
equation  are  known  except  K0     K  can  be  determined  by  setting  the 
equation  equal  to  a  very  small  number  and  solving  for  K„     The 
error  in     P  (Z-n)     will  then  be  less  than  that  small  number  times 

p     e         <^l=  o 


that  the  value  of  K  has  been  approximated ,  the  probability 

distribution  of  demand  during  lead  time  can  be  found  as 

K 


In  order  to  simplify  the  calculation  of  this  number  manually  or  in 
a  computer 8   the  technique  of  "nesting"  will  be  shown,, 
K 
I    f*"1)     ak(k+Dn  =  1+Q  a2n+f+1)   a2  3°  +   ...  +    (r+^1)aK(K+l)n 


k=0 


«  1.  +  »  2n+  £%tH     a2  3n  +  ...  ♦  r(rH)...(r+K-l)  &K(K+l)i 

2 1  I\! 


•('*-(«¥•  r'M-^-fT)"  ))•••)) 


This  can  then  be  written  in  the  compact  form. 

Let 

r+K=1        /  K+1   \n 
A, 


r+K=1        /  K+1   V 
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and 


30  =  1   +  AQ 


B1  s  1   +  A1   B0 


B.  =  1   +  A.  B,   „ 


where 


r+K-1-1       /K-j+1\ 

A  ,  =  -— ^'-—  a  [      J       ]    and  B   .   =  1 

\   K-.1    )  ~1 


J  K-j  V   K-j 

then 


P   |z  =  n]=[l+*K_1BK_a].-,kpri|. 


n 
nl 

As  another  example 
Example  2 
Let 

D  be  a  negative  binomial  random  variable 
L  be  a  translated  Poisson  random  variable 
Then 

f(s)  s[— M  from  eq,    U.3) 

\1=qs/ 


and 


\  A(s=1)  ,       . 

;)  =  s  e  from  eq0    (4.. 2) 


Then 


h(s)  =  f(s)  g  [f(s)]   =  f(i)  eXCf(s)-1] 


2U 


Let  us  again  find  the  Maclaurin  expansion 


hw  =  feMfef  -1H 


CO 


A      r  (k+1)   (.^a\-r(k-¥l) 

£6  T?  p  (1=qs) 


CO  CO 

2,    A       r(kH)    A    /r(k+l)+n-l\    /     a)n 
^  ~k7  P  n=0   V        n  j    (q  S) 


V      n     -*      n  V  /r(k+l)+n-l\     A       r  (k+1) 
2L    s     e         q      ^  7:  P 

n^O  k=0  V         n  '      k! 


and  od 

(APr) 


p    z 


1        ->      n     r  Y    /r(k+D+n-l\ 
-  n J  -  e         q     p     Z. .1  n  j 


k=0    *  k! 


To  find  a  bound  for  the  remainder  of  the  series 
Let 

K  K 


(APr)k 


k-0  k=0  \         n  / 


and 

CO  00 


Z,      %  ^    %   Z,  "T 


%  "  ,  *->      "k  "      "K  f—t      u 

k=K+1  1=0       K 


Then 


l(ktlj±n.lll,      ak  where  a  =    ^ 

k       (r(k+l)-1j   I    n! 
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and 


UK-:  (r  K+rl  +  n~1 ) !      (rK  +  r-1)    !        i 


UK  (p  K  +  pl=1 ) !    (r  K  +  p+n-1 ) ! 


-     (rK+rl+r+c-1 )    (rK+rl+n-2) . . .  (rK+r+n)       1 
(rK+rl+r~l)    (rK+rl+r-2)    60C   (rK+r) 


4  v  +  )  (i  + 

V  p(K+1+1)-1    /    V 


p(K+1+1)-2 


I         V  r(K+lW 


Jl 


^  v      p(k+i); 


Thus 


UK+1 


* 


r  -il 


fi+^5Jo)  a] 


let 


*-  [(1+^r)rtl 


We  then  have  to  choose  minimum  K  such  that  Y<  1  so  the  resulting 
geometric  series  can  be  summed,  i.eOJ 


*=[(i+^n<1 
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-If) 


if 

1  +  - 

and 


,(D 


v\n    &     


Then        X  <  1 
Now 

CO  CD 

y  °k+i  .    y 

K  K   1=0  K  K    1=0 


=  uK  = 


i-r 


and  to  solve  for  K  we  have 


-*-■-*         I     y\  _    &  (W)-HM]   1     prK  «V  Pf 


e       q     p     uK    - 


[r   (K+1)   -l]!n!    [l-  p*   (1+  %-  )r] 


Upon  setting  this  equation  equal  to  an  allowable  error,  K  can  be 
determined o 

The  estimate  of  the  remainder  is  the  value  of  the  allowable 
error  from  which  K  was  determined. 
Then,  given  K 


K 


_]T    Mk+D+n-l\      Upr)k 
"k=0\         n  /         k! 
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The  "nesting"   technique  is  helpful  in  computing  the  probabil- 
ities  of  demand  during  lead  timeu 


«.     _  ^"°^  ^r  _n  « 
n     =  e       p    q     b.% 


,  r  >1 


and 


AR  = 


.(r+n-lll  +  _(2rtn-l)j,.(Apr)  +  (3r+n-l) !  (Apr)2  + 
(p~l)ln!  (2r-l)!n!  (3r-l)l    n! 


-  (rfJ^lJl  U  +  (2r+n°l)oc.(r+n)    (Apr)      (  1  + 
(p-l)!nl   V  (2r-l)      ...    (r) 


U  +     (rk+n-1)    ...    (rk+n-r)    (Apr) 
V  (rk-1)   ...     (rk-r)k 


-f-      o  o  o 


..A  +     (rK+n-1)...    (rK  +n-r)    (A  pr)     \\       tt 
\  (rK-1)    ...      (rK-p)  //        // 


Thus,  if 


Bj      L   +     r(K-j)-lJ  L        p(K-j)-2j  L         r(K-j)-p_|AP 


and 


Cr>  =  1  .  then 


C   = 


'j+1 


1  +  B0  C0 


1   +  Vi 
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and 


1  +  BK-2  CK-2 


Now  that  the  distribution  of  demand  during  lead  time  and  methods 
for  calculating  the  probabilities  have  been  discussed,  it  is  a  simple 
matter  to  find  the  number  of  items  of  stock  that  are  required  as  the 
reorder  point.  The  reorder  point  must  be  the  number  of  units  of  stock 
that  will  last  until  new  replacement  stocks  arrive  with  the  probabil- 
ity of  a  stockout  during  the  lead  time  small. 

Let  a  =  Acceptable  probability  of  a  stockout  during  lead  time 


Then 

P  |Z>Rj<a 
must  be  solved  for  the  value  of  R. 


[z>h] 


i-p[z>rj  =  p   z  =  q+p[z  =  i]+...  +  p[z  =  R] 

o " .  R  is  the  smallest  number  such  that 
R 

2_.P|_Z=n|^1-a 
n=0 
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5.       Negative  Binomial  Distribution  of  Demand  During  Lead  Time 

It  would  be  very  helpful  if  the  distribution  function  of  de- 
mand during  lead  time  could  be  found  for  any  distribution  functions 
of  demand  and  lead  time  in  a  familiar  form,  but  this  is  not  possible, 
It  can  be  done  in  some  special  cases.     One  known  case  is  given  below 
as  an  example.,     Consider  deriving  the  mass  function  of  demand  during 
lead  time  where? 

(1)  The  demand  is  random  and  the  number  of  demands  d 
during  a  fixed  interval  of  time  t  has  a  Poisson 
distribution* 

=at       fo+^d 
p(djat)  *  *  vaTV  ,d  =  0,1,2,... 

where  a  is  the  average  number  of  demands  per  unit  time. 

(2)  The  lead  time,  t,  has  a  gamma  type  probability  density 
function 


~ut       /    .\k-1 
f(t,u,k)  m&  l™  ,  t>0 


k  k 

with  mean  —*-.  and  variance  — — 

u  n* 


Then  the  probability  of  exactly  z  demands  during  the  lead  time  is 


JCO 


^-at       /   iNz  .__-ut       /..iNk-1 

p(z$k,u,a)  = 


(at)  ue  (ut) 


z!  (k-1) 
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.  f      z+k-1     -(a+u)t 

az  uk  t  e  dt 

!(k-l)!        o 


Letting    x  =  (a+u)   t,  then 


z       k 
a       u 


p(z^k,u,a)  =  7 r  z+k       ;       T      /      xz+k~1   e~x     dx 

*     '    *    s  (a+u)z+K  z!(k-l)!   ^o 


az  uk  (k+e-1)! 


(a+u)Z+k       z!(k-l)! 


k+z~1 

z 


(-fry 

\a+u/     \a+u/ 


the  last  step  being  valid  if  k  it   an  integer ,     Thus  the  demand  during 
lead  time  has  a  negative  binomial  distribution  with 

ak 
mean  =  =— 


and 


variance 


ak  /         a\ 

x-       1    +- 

U      \  U/ 
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6.   Summary  and  Conclusions 

The  distribution  of  demand  during  lead  time  has  been  discussed 
when  discrete  random  variables  describe  the  demand  and  lead  timee 
It  is  unfortunate  that  the  distributions  do  not  appear  as  well 
known  distributions,,  The  use  of  probability  generating  functions 
has  been  shown  as  a  method  by  which  the  distribution  can  be  found 
and  used.  Calculations  of  the  probabilities  are  cumbersome  since 
in  the  examples  shown  each  of  coefficients  of  the  s  in  the  prob- 
ability generating  function  involved  an  infinite  sum,,  The  amount 
of  computer  time  required  for  the  calculations  has  not  been  investi- 
gated, but  since  the  calculations  are  straight  forward  and  the  number 
of  terms  required  before  the  infinite  sum  converges  to  an  approximate 
value  are  relatively  small,  the  use  of  the  procedure  outlined  in 
Section  U   should  be  useful  and  not  too  cumbersome  to  preclude  their 
use. 

The  use  of  the  generating  function  method  requires  data  for  the 
estimation  of  the  parameters  in  the  demand  and  lead  time  distribu- 
tions. This  data  is  readily  available  in  the  historical  files  of 
an  inventory  system.  Other  methods  of  estimating  the  demand  during 
lead  time  requires  collection  of  data  directly  from  the  number  of 
demands  during  a  lead  time.  This  data  is  not  usually  collected  and 
even  if  it  were  available  sufficient  data  to  obtain  a  good  statis- 
tical sample  would  require  a  very  long  time. 

Further  investigation  of  distributions  calculated  by  the  methods 
in  Section  4-  may  show  that  the  probability  distributions  could  be  fit 
to  some  known  distributions  and  the  resulting  distribution  used  to 
approximate  the  demand  during  lead  time. 
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This  Investigation  has  been  made  primarily  to  consider  a  method 
of  determining  the  demand  during  lead  time  distributions  for  low 
demand  items ,  Low  demand  items  are  particularly  troublesome  in  an 
inventory  system0  Military  spare  parts  inventories  have  large 
numbers  of  very  low  demand  i terns 0  The  low  demand  items  are  also 
characterized  by  high  variances  [6], 

The  use  of  the  Poisson  and  negative  binomial  distributions 
have  been  used  extensively  in  this  paper  since  these  distributions 
have  been  found  to  fit  low  and  medium  demand  items .  The  Poisson 
distribution  is  generally  used  for  items  with  a  demand  of  two  or 
less  items  per  year.  The  negative  binomial  distribution  has  been 
used  for  items  with  demands  ranging  from  two  to  one  hundred  per 
year  [7]0  In  the  case  where  the  Poisson  demand  distribution  is 
applicable,  the  generating  function  method  for  establishing  a  re- 
order point  can  be  readily  applied.  Since  the  reorder  point  will  be 
relatively  low,  the  number  of  probabilities  that  must  be  calculated 
before  the  sum  is  equal  to  or  greater  than  the  service  level  desired 
will  be  small ,  It  has  not  been  discussed  in  this  paper,  but  when 
reorder  points  have  been  previously  established,  the  probability 
generating  function  for  demand  during  lead  time  can  be  used  to  find 
the  tail  of  the  distribution.  With  a  reorder  point  previously  estab- 
lished, a  good  estimate  of  the  value  of  reorder  point  is  available 
and  with  several  trials  on  either  side  of  the  old  reorder  point, 
the  tail  of  the  distribution  could  be  used  to  determine  the  new  order 
point.  In  this  manner,  the  probabilities  would  not  have  to  be  summed 
until  the  sum  was  greater  than  the  desired  probability.  Feller  [2] 
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gives  the  following  relations  between  the  probability  generating 
function  and  the  generating  function  of  the  tail  of  the  distri- 
bution. 

Let  Q  (s)  =  the  generating  function  of  the  tail 


Then 

1-  h(s) 


Q  (s)  = 

1-s 


The  calculation  of  the  tail  in  this  paper  has  not  been  investigated 
since  it  gets  more  complex  than  the  probability  generating  function. 
However,  in  special  cases  and  in  a  particular  situation  as  described 
above,  the  tail  of  the  distribution  may  be  the  easiest  method  for 
setting  a  reorder  point. 

The  special  case  of  Poisson  type  demand  and  gamma  distributed 
lead  times  was  discussed  because  this  combination  of  demand  and  lead 
time  distributions  lead  to  a  negative  binomial  distribution  of  de- 
mand during  lead  time.  Further  investigation  of  this  model  with 
actual  data  may  be  very  fruitful.  The  Poisson  distribution  is  a 
good  distribution  for  low  demand  items  and  the  gamma  distribution 
is  a  very  versatile  two  parameter  distribution  that  can  be  used  to 
approximate  a  wide  range  of  lead  time  patterns  [5].  Consequently, 
any  demand  and  lead  time  process  that  may  possibly  be  fit  to  this 
model  should  be  considered,  especially  since  the  probability  law  of 
demand  during  lead  time  may  be  readily  computed. 

The  writer  wishes  to  thank  Professor  Robert  R.  Read  of  the 
U.S.  Naval  Postgraduate  School  for  his  suggestions,  assistance,  en- 
couragement and  cooperation  in  the  preparation  of  this  paper. 
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APPENDIX  I 
MEAN  AND  VARIANCE  FROM  PROBABILITY  GENERATING  FUNCTIONS 
The  mean  and  variance  of  a  discrete  random  variable  may  be 

calculated  from  the  probability  generating  function.  It  may  be 

shown  that  (see  e.g.,  Feller  [2]) 

E(X)  =  f'(l) 

Var(X)  =  f"(l)  +  f'(l)  -  [f'(l)]2 


As  an  example ,  the  mean  and  variance  of  demand  during  lead 
time  when  lead  time  is  a  translated  Poisson  random  variable 

x-1 


p(x|A)  =  e' 


■a  _A 


9   X  —  1  97 , 


(x-1)! 

with  mean  =  (A  +1 )  and  variance  =  A 
and  demand  has  the  negative  binomial  distribution 

p(x;rsp)  =  ^"  *  j  pr  qx      ,  x=0,1,2,°" 

with  mean  =  -—-*  and  variance  =  — j 
P       .        P 

then  the  random  variable  Z  has  a  probability  generating  function 
h(s)  =  | 


1-qs 


Kte-r  •']) 


taking  the  first  derivative  with  respect  to  s  we  have 

T+1  \2 

-qs, 


h.(s)  =  ^(-^.p^^fy-l    +  ^  A(^fr+1ex{A[[^| 
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and 


h!(l)  =  —2-  +  -— 


a-d+A) 


the  second  derivative  evaluated  at  one  is 

h««(l)  =  -ij  f[rq(l+A  )]2  +  r2q2^  +  rq2(l+A)| 

then 

Var(Z)  =  h"(l)  +  h'd)   -  [h»(D] 


"   (1+A)  +  (? J" 


P2 


These  results  can  be  verified  against  the  general  formulas  for 
the  mean  and  variance  given  in  Section  2. 
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